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ABSTRACT 



Aims. We search for low-mass companions in the innermost region (< 300 mas, i.e., 6 AU) of the j} Pic planetary system. 

Methods. We obtained interferometric closure phase measurements in the K-band with the VLTI/AMBER instrument used in its 

medium spectral resolution mode. Fringe stabilization was provided by the FINITO fringe tracker. 

Results. In a search region of between 2 and 60 mas in radius, our observations exclude at 3<r significance the presence of companions 
with K-band contrasts greater than 5 x 10~ 3 for 90% of the possible positions in the search zone (i.e., 90% completeness). The median 
Icr error bar in the contrast of potential companions within our search region is 1.2 x 10" 3 . The best fit to our data set using a binary 
model is found for a faint companion located at about 14.4 mas from ji Pic, which has a contrast of 1.8 x 10~ 3 ± 1.1 x 10~ 3 (a 
result consistent with the absence of companions). For angular separations larger than 60 mas, both time smearing and field-of-view 
limitations reduce the sensitivity. 

Conclusions. We can exclude the presence of brown dwarfs with masses higher than 29 Mj up (resp. 47 Mj up ) at a 50% (resp. 90%) 
completeness level within the first few AUs around /? Pic. Interferometric closure phases offer a promising way to directly image 
low-mass companions in the close environment of nearby young stars. 

Key words. Stars: individual: /? Pic - Planetary systems - Planets and satellites: detection - Techniques: interferometric 



1. Introduction 



The young (~ 12 Mvr. lZuckerman et al.ll200lT) . nearby (19.3 pc), 
and bright (K = 3.5) A5V-type star ft Pictoris (HD 39060) 
is surrounded by one of the most famous extrasolar planetary 
systems, consisting of a rece ntly detected planetary companion 
dLagrange et al . 2009a, 2010) i nside a n optically thin debris disk 
seen edge-on (Smi th & Terrildfl984l) . which has been resolved 
at various wavelengths. Several asymmetries have b een identi- 
fied in the debris disk, including a warp at ~ 50 AU dHeap et alJ 
2000) that is now understood to be the result of the dynamical 
influence of a massive body (a fe w Jupiter masses) on a n eccen- 
tric orbit around the central star (Freis tetter et alJl2007|) . The 9- 
Mj up companion discovered by Lagrange et al. (2009a) may be 
the cause of this warp. We note that the planetary nature of thi s 
companion was not easy to ascertain (ILagrange et alJ 120 09b), 
because the companion was located at a projected distance 
smaller than the inner working angle of VLT/NACO (335 mas 
for 9 Mj up ) between the discovery observations in 2003 and the 
confirmation observations in late 2009. 

Long-baseline optical interferometry is a promising tech- 
nique to search for faint companions at angular separations 
smaller than the diffraction limit of a single aperture. In partic- 
ular, closure phase measurements on a closed triangle of base- 
lines are very sensitive to asymmetries in the brightness distri- 



* Based on observations collected at the ESO La Silla Paranal 
Observatory under program IDs 084.C-0566 and 384.C-0806. 
** FNRS Postdoctoral Researcher, email: absil@astro.ulg.ac.be 



bution of the source, and can be used to detect faint companions. 
Closure phases have the added advantage of being insensitive to 
telescope-specific phase errors (unlike visibilities and phases), 
including a tmospheric tur bulence effects (for a review of closure 
phases, see lMonnieij|2003l) . 

The interferometric detection of extrasolar planets (hot 
Jupiters in particular) has already been attempted by a few 
groups, using in particular precision closure phase measure- 
ments. Despite the exquisite accuracy that has already been 
reached (e.g., 0!1 stability in the CHARA/MIRC closure phases, 
Zha o et al. l T2008h . no extrasolar planet has yet been detected. 
Differential phase techniques have not been more success- 
ful because of atmospheric and instru mental limitations (e.g., 
iMillour et all 120081: [Matter et alJ I201CJ) . Higher sensitivities to 
faint companions can be reached when closure pha ses are ob- 
tained on fully resolved stellar photospheres (e.g.. lLacour et al.l 
2008; iDuvert et alJl20Toh . but this is unfortunately not the case 
for most main-sequence stars with currently available interfer- 
ometric baselines. In this Letter, we perform a deep interfero- 
metric search for faint companions at short angular distances 
from the unresol ved young main seq uence star ft Pic (6>ld = 
0.85 + 0.12 mas. iDi Folco et alJ. 12004). based on closure phase 
measurements with VLTI/AMBER. 

2. Observations and data reduction 

Observations of ft Pic were performed on four different nights 
from 2010 January 24 to 28 with the AMBER instrument, used 
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in its medium resolution mode (R =s 1500 from 1.93 to 2.27 /im) 
with three 1.8-m Auxiliary Telescopes of the VLTI arranged on 
the A0-G1-K0 triangle (ground baselines from 90 m to 128 m). 
Fringe tracking was provided by the FINITO facility during all 
our observations. A total of 12 Observing Blocks (OBs) were 
completed for B Pic, with Detector Integration Times (DIT) 
ranging from 0.5 to 1 sec depending on the atmospheric condi- 
tions. To reduce time smearing in our data set, long OBs were 
divided into shorter sequence of 15 min or less, giving a to- 
tal of 26 OBs. Our observations of B Pic were interleaved with 
observations of HD 39640, a G8III calibration star of magni- 
tude K — 3.0 and a ngular diameter 6>ld = 1.235 + 0.015 mas 
( Mera ncTet alJ l2005). located only 1°2 from/3 Pic on the sky. 

Raw data were reduced using the amdlib v3.0 package 
riTatulli et alJl2007l: IChelli et al.ll2009h . using all recorded detec- 
tor frames (no fringe selection). The B Pic closure phases were 
calibrated by a simple subtraction of the average closure phase of 
all calibrator measurements obtained during the same night with 
the same instrumental set-up. This procedure implicitly assumes 
that the closure phase transfer function does not vary signifi- 
cantly during several hours of observations, a hypothesis that we 
verified during our four observing nights. The error bar related 
to the calibration process is generally estimated by measuring 
the standard deviation in the calibrator closure phases during the 
night. In our case, too few calibrator measurements are available 
to provide a reliable error bar, so we decided to use the standard 
deviation in all measurements (science and calibrator) instead. It 
must be noted that, if the scientific target were to have a com- 
panion producing a significant closure phase signal, this would 
lead to an overestimated calibration error bar and could possi- 
bly hide the companion (we address this in the last paragraph 
of Sect. 13. 11 1. The evaluation of the calibration error bar is illus- 
trated in Fig.[T]for the night of January 24, which is actually the 
poorest night in our data set in terms of transfer function sta- 
bility. In this plot, the data were binned into 6 spectral channels 
to reduce their statistical errors and isolate the effect of trans- 
fer function instabilities. The short end of the K band (top panel 
of Fig. Q]) displays significant instabilities in its closure phases, 
most probably caused by strong and variable spectral lines in 
the Earth atmosphere. These instabilities disappear beyond about 
2.0 yum, so we restrict our spectral range to the 2.00-2.26 /vm re- 
gion. In this spectral range, we find that the RMS calibration er- 
rors range between (f.20 and (f.31 depending on the night. These 
error bars will be added quadratically for our whole data set, 
night by night (after binning individual spectral channels, where 
required). 

3. Searching for faint off-axis companions 

3.1. Data analysis 

The field-of-view (FOV) of the AMBER instrument is limited 
by the use of single-mode fibers. We estimated the off-axis trans- 
mission of a point-like object in the presence of atmospheric tur- 
bulence corrected for tip-tilt fluctuations by the STRAP system 
on the ATs. Our simulations produced a Gaussian transmission 
with a full width at half maximum of 420 mas for the median 
seeing of our observations (~ 078). For angular distances larger 
than 210 mas, the transmission quickly drops to low values and 
we limit our analysis to separations smaller than about 300 mas, 
where the transmission drops by a factor of 4. This corresponds 
to a linear distance of ~ 6 AU from B Pic. We note however 
that, because of the Earth's rotation, closure phases are variable 
in time. In order not to smear the closure phase signal too much 
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Fig. 1. Fluctuation of the closure phase as a function of time for 
/3 Pic (black empty diamonds) and its calibrator (blue filled di- 
amonds) during the night of January 24. Data have been binned 
into 6 spectral channels to reduce statistical errors. The solid 
lines represent the mean of the CP transfer functions and the 
dashed lines their standard deviations. 



during our ~15 min OBs, we had to restrict our search region to 
a FOV radius of about 50 mas. A finer time sampling would be 
needed to explore the full 300 mas FOV at the highest sensitivity. 

The first step in our data analysis is to adapt the spectral res- 
olution to the explored FOV, by ensuring that the variations in 
the closure phase as a function of wavelength, created by a po- 
tential companion, are sampled with at least four data points per 
period. The periodicity in the closure phase signal as a function 
of wavelength is given roughly by = A 2 l{B6 - A), where B is 
the mean interferometric baseline length. For B 100 m in our 
case, and a maximum angular separation of 50 mas, this implies 
that a; 0.2 fim. The appropriate spectral bin size to cover a 
50 mas FOV is thus AA = P A /4 = 50 nm, so that a total of 5 
spectral channels are needed across the 2.00 - 2.26 //m spectral 
range. We used a 5-sigma clipping method to remove outliers 
when binning our data into the synthetic spectral channels, and 
statistical error bars are estimated to be the standard deviation in 
the individual data points of each bin. 

The next step is to search the entire FOV for possible com- 
panions. To do this, we used th e photospheric model of B Pic pro- 
posed bv lDi Folco et al . (2004), to which we added a secondary 
point-like companion at various locations within the FOV. The 
X 1 distance between the data and all our binary models was then 
computed, and for each position we selected the companion con- 
trast that produces the lowest^ 2 . The resulting^ 2 map was then 
converted into a probability map (Pig. [2j using the^ 2 probability 
distribution function, taking into account the number of degrees 
of freedom in our x 1 distribution. In the present case, there are 
26 x 5 independent data points and 1 parameter to fit (the com- 
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Fig. 2. Probability of a binary model to reproduce our data set, for various positions of the secondary companion across a FOV 
50 mas in radius. At all places in the FOV, the flux of the companion has been tuned to minimize the x 2 -- The zones in white have 
their minimum xl for a com panion of zero flux. T he position of /3 Pic is represented by a star, and the orientation of the circumstellar 
disk midplane (PA = 29°5. lBoccaletti e t al. 2009) by a black dotted line. The u, v coverage of our data set is represented in the upper 
right inset, including wavelength dependence (using colors from blue to red), with the axes graduated in cycles per arcsec. 



panion contrast for each binary position), hence 129 degrees of 
freedom. In the absence of a companion, the reduced chi square 
(X 2 r ) amounts to 1.01, which corresponds to a probability of 45% 
for our photospheric model alone to reproduce the data set. 

Our probability map shows a local maximum of 86% (i.e., 
X 2 r = 0.87) at a position (ARA, ADec) = (3.9, -13.9) mas, with 
a best-fit contrast of 1 .8 x 10~ 3 . The quality of our best-fit model, 
illustrated in Fig. [3] confirms that this model closely reproduces 
our data set. Several other companion positions within our FOV 
would provide (almost) equally good fits. We note however that 
the pure photospheric model also reproduces the data quite well, 
so that there is no real evidence of a companion. We now verify 
whether this could be caused by an overestimation of the calibra- 
tion error bars, as discussed in Sect. [2] We artificially decrease 
our estimated calibration error bars until we achieve y 2 r — 1 for 
our best-fit model (which is reached when calibration error bars 
are multiplied by 0.92). The significance of this possible detec- 
tion can then be estimated by converting the probability of the 
null hypothesis (i.e., no companion), which is now only 9%, into 
an equivalent number of standard deviations for the possible de- 
tection of a companion, using standard relationships for normal 
distributions. We infer a 1 .7<x significance for our detection, giv- 
ing a companion contrast of 1.8 x 10~ 3 ± 1.1 x 10~ 3 . This low 
significance confirms that it cannot be considered as a real de- 
tection. 



3.2. Sensitivity limits 

Sensitivity limits can be derived from our xl ma P (° r equiva- 
lently, probability map), by searching at each point of the FOV 
for the companion contrast that would produce a xl larger than 
a pre-defined threshold. We choose a 3<x criterion to define 
our sensitivity limit (i.e., probability of less than 0.27% for the 
model to reproduce the data). With 129 degrees of freedom, this 
corresponds to a threshold X 2 -u m - 1-38. Once the 3cr upper limit 
to the detectable companion contrast is computed at each point 
of the FOV, one can define global sensitivity limits by building 
and inspecting the histogram of the 3<x detection levels across 
the considered FOV. For instance, for the 50 mas FOV used in 
Fig. [2] the contrast upper limit is < 3.6 x 10~ 3 for 50% of all 
positions in the FOV (50% completeness) and < 4.8 x 10~ 3 for 
90% completeness. This confirms that the typical lcr error bar in 
the contrast of a detected companion would be about 1.2 x 10~ 3 , 
as estimated in the previous paragraph. 

These sensitivity limits were confirmed by a double-blind 
test, where synthetic companions of various contrasts were in- 
troduced into our raw data set. In these blind tests, we were able 
to retrieve the companions with a contrast of 3.0 x 10~ 3 in about 
50% of the cases (although with a formal significance gener- 
ally between 2 and 3cr), and the companions with a contrast of 
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Fig. 3. Representation of the whole data set (diamonds with error 
bars, binned into 5 spectral channels) and of the best-fit model 
(red curves) found on the 50 mas radius FOV. This model corre- 
sponds to a companion with a contrast of 1.8 x 10~ 3 located at 
cartesian coordinates ARA = 3.9 mas and ADec = -13.9 mas 
relative to the central star. 



5.0 x 10 -3 in all cases. These tests confirm the validity of our 
contrast upper limits based on the^ analysis. 

Sensitivity limits can also be computed as a function of an- 
gular separation, by building x 2 r maps on annular fields-of-view 
of increasing size. This is illustrated in Fig. [4] where annuli 
10% in relative width have been used. This figure shows that 
VLTI/ AMBER reaches its optimum sensitivity in the 2-60 mas 
region, where a median contrast of 3.5 x 10~ 3 ( 5.0 x 10~ 3 ) can 
be reached at a 50% (90%) completeness level using the ATs. 
Beyond 60 mas, the effect of time smearing on the closure phase 
signal of potential companions becomes significant, reducing the 
sensitivity. For larger separations (> 200 mas), the sensitivity 
degrades more rapidly because of the decreasing off-axis trans- 
mission of the single-mode fibers. The inner working angle of 
our interferometric study is about 1 mas, where the sensitivity 
drops to a few percent in contrast (e.g., 5 x 10~ 2 at 90% com- 
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Fig. 4. Sensitivity curves showing the 3cr upper limit to the con- 
trast of off-axis companions as a function of the angular sepa- 
ration for 50% and 90% completeness, computed across annu- 
lar fields-of-view with 10% relative widt h. Equivalent masse s 
were computed using the COND model of Baraf fe et al.l (|2003), 
for an age of 12 Myr. The companion discove ry zones of ra- 
dial velocity measurements dGalland et al.l l2006. left of the blue 
dash-dotted line) and of AO-assisted coronagraphic imaging 
dBoccaletti et al.ll2009l right of the green dash-dotted line) are 
shown for comparison. 



pleteness) because of the limited angular resolution provided by 
our baselines. 



4. Discussion 

To compare our sensitivity limits with other studies, it is useful 
to express them in terms of companion masses. The first step 
is to convert contrasts into absolute magnitudes, taking into ac- 
count the K magnitude and distance of /3 Pic. Absolute magni- 
tudes are then converted into masse s using the COND evolution- 
ary models of Bar affe et al.l d2003l) . assuming an age of 12 Myr. 
The result is illustrated in Fig. [4] where mass upper limits are 
represented by dotted lines. In the 2-60 mas region, the median 
upper limit to companion masses is 29 Mj llp (47 Mj up ) for 50% 
(90%) completeness. 

These upper limits should be compared with the main planet 
search methods around nearby main-sequence stars: radial ve- 
locity (RV) monitoring and direct (single-pupil) imaging. State- 
of-the-art RV meas urements obtained with HARPS on /3 Pic 
dGalland et alj|2006l) have reached an accuracy of 180 m s _1 (af- 
ter correction for its pulsations), which provides a typical sen- 
sitivity of 10 Mj up for a semi-major axis a = 1 AU. The sen- 
sitivity then scales as a~ 2 . The largest semi-major axis that can 
be reached depends on the timescale of the RV monitoring: to 
cover most of our interferometric search region (up to 6 AU in 
semi-major axis), j3 Pic should be surveyed for about 5 years. 
For this time coverage, the only companions that could be de- 
tected by our interferometric search and not by the RV mon- 
itoring would be those located at orbital distances larger than 
6 AU, which would by chance be at projected angular separa- 
tions smaller than 300 mas at the time of our observations. K- 
band direct imaging observations using a Four Quadrant Phase 
Mask (FQPM) coronagraph on VLT/NACO have yielded an up- 
per limit of 2.5 x 10~ 4 to the con trast of off-axis compa nions 
at projected distances > 335 mas (Lagrange et al. 2009b), and 
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can reach a cont rast of 1.3 x 10~ 2 at t he FQPM inner working 
angle of 70 mas dBoccaletti et al.l l2009). In practice, the NACO- 
FQPM sensitivity becomes superior to that of AMBER for an- 
gular distances larger than ~ 100 mas, so the two techniques can 
be considered complementary. 

Although our detection limits are promising and would al- 
low brown dwarfs to be detected around bright young stars such 
as j6 Pic, they are insufficient to detect the planetary companion 
discovered by Lagrange et al.l d20b9a). which has an estimated 
K-band contrast of 2.5 x 10~ 4 . The accuracy of our measure- 
ments needs to be improved by at least a factor of 1 to reach this 
level of performance. The same improvement in measurement 
accuracy would be required to bring our sensitivity down to the 
realm of hot- Jupiter type planets. It is expected that the next gen- 
eration of interferometric instruments at the VLTI (PIONIER, 
GRAVITY, and eventually VSI) could provide the necessary im- 
provements in instrumental stability and sensitivity to achieve 
this performance. 
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